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New data from parity-violating experiments on the deuteron now allow isolation of 
\Q • the strange-quark contribution to the nucleon magnetic moment, G%^{0), without the 
^ ■ uncertainty surrounding the anapole moment of the nucleon. Still, best estimates place 
^m(O) > 0- It is illustrated how this experimental result challenges the very cornerstone 
^ i of the constituent quark model. The chiral physics giving rise to G%i{0) ~ is illustrated. 

in 

O ■ 1. INTRODUCTION 

^ ■ Previous experimental reports [ of the strange-quark contribution to the nucleon 
Q . magnetic moment, G^(0), have involved the poorly understood nucleon anapole moment. 

However, combining the original proton data with new data from parity violation exper- 
iments on the deuteron now allow an independent isolation of Glj{0). The preliminary 
CJ I analysis [ H indicates 6*1/(0) is smaller than the earlier result [[^ of 
^ '. 

?h; ^^^(0.1 GeV^) = +0.61 ±0.17 ±0.21/ijv, (1) 



X 



but positive values for 6*^/(0) are still favoured. 



;_j ■ In this report, we elucidate the impact of this experimental result on our understanding 
of the concept of a constituent quark. In particular, the constituent quark approach places 
properties such as the magnetic moment contribution of a quark to a baryon intrinsic to 
the constituent quark itself and therefore not dependent upon the environment in which 
the quark resides. 

The extraction of G^(0) is based on the approximation of charge symmetry in the 
nucleon moments; i.e. equal u and d current quark masses and negligible electromagnetic 
corrections. This approximation is expected to be valid at the 1% level [ Q and is certainly 
sufficient for the analysis of (7^/(0). 

An examination of the symmetries manifest in the QCD path integral for current matrix 
elements reveals various relationships among the quark sector contributions [ Q . Current 
matrix elements of hadrons, such as the magnetic moment of the proton, are extracted 
from the three-point function, a time-ordered product of three operators[ Generally, an 
operator exciting the hadron of interest from the QCD vacuum is followed by the current of 
interest, which in turn is followed by an operator annihilating the hadron back to the QCD 
vacuum. In calculating the three point function, one encounters two topologically distinct 
ways of performing the electromagnetic current insertion. Figure |l| displays skeleton 

*Presented at Few Body 2000, XVIth lUPAP International Conference on Few-Body Problems in Physics, 
National Taiwan University, Taipei, Taiwan, March 6-10, 2000. 



2 




g g 



Figure 1. Diagrams illustrating the two topologically different insertions of the current 
within the framework of lattice QCD. These skeleton diagrams for the connected (left) 
and disconnected (right) current insertions may be dressed by an arbitrary number of 
gluons and quark loops. 

diagrams for these two possible insertions (with Euclidean time increasing to the right). 
These diagrams may be dressed with an arbitrary number of gluons and quark loops. 
The left-hand diagram illustrates the connected insertion of the current to one of the 
"valence" quarks of the baryon. In the right-hand diagram the external field produces 
a. qq pair which in turn interacts with the valence quarks of the baryon via gluons. It is 
important to realize that within the lattice QCD calculation of the loop diagram on the 
right in Fig. |l] there is no antisymmetrization (Pauli blocking) of the quark in the loop 
with the valence quarks. For this reason, in general only the sum of the two processes in 
Fig. |I| is physical. 

2. CHARGE SYMMETRY EQUALITIES 

Under the assumption of charge symmetry, the three-point correlation functions for 
octet baryons leads to the following equalities for electromagnetic current matrix elements 

p= euUp + eddp + On , n = e^Up + e^dp + On , 

= Cs 5= + e„ Mho + 0= , S~ = 5= + Cd u^o + 0= . (2) 

Here, O denotes the contributions from the quark-loop sector - shown on the right-hand 
side of Fig. |l]. The baryon label represents the magnetic moment. Subscripts allow for 
environment sensitivity implicit in the three-point function For example, the three- 
point function for E"*" is the same as for the proton, but with d replaced by s. Hence, the 
M-quark propagators in the S"*" are multiplied by an s-quark propagator, whereas in the 

^Note that the term "valence" used here differs from that commonly used in the framework of deep- 
inelastic scattering. Here "valence" simply describes the quark whose quark flow line runs continuously 
from — > a;2. These lines can flow backwards as well as forwards in time and therefore have a sea 
contribution associated with them [ || . 
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proton the ti-quark propagators are multiphed by a (i-quark propagator. The different 
mass of the neighboring quark gives rise to an environment sensitivity in the w-quark 
contributions to observables, which means that the naive expectations of the constituent 
quark model Up/uY,+ = Un/usp = 1 may not be satisfied [ ^, ^ |, ^, [1^, |11|]. This 
observation should be contrasted with the common assumption that the quark magnetic 
moment is an intrinsic-quark property which is independent of the quark's environment. 

Focusing now on the nucleon, we note that for magnetic properties, contains sea- 
quark-loop contributions from primarily u, d, and s quarks. In the SU(3)-fiavor limit 
{rriu = rrid = rris) the charges add to zero and hence the sum vanishes. However, the 
heavier strange quark mass allows for a result which is non-zero. By definition 
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where R = 



(3) 
(4) 



and the leading superscript, £, reminds the reader that the contributions are loop con- 
tributions. Note that, in deriving Eq.(|D, we have set ^GJ{,j = ^Gfj, corresponding to 
mu = rud [ §]. In the constituent quark model ^R^ = md/rris — 0.65. 

With no more than a little accounting, the strange-quark loop contributions to the 
nucleon magnetic moment, G\.j may be isolated from (j^) and in the following two 
phenomenologically useful forms. 
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(5) 



(6) 



As we explained above, under the assumption that quark magnetic moments are not 
environment dependent, these ratios (i.e. Up/uY,+ and m„/m=o) are taken to be one in 
many quark models. Incorporating the experimentally measured baryon moments leads 
to: 



M 



G 
and 



R 



l-'R 



3.673 



ns4 



(3.618) 



'R 



-1.033 - 



-0.599) 



(7) 



(8) 



where all moments are expressed in nuclear magnetons (/iat). (Note that the measured 
magnetic moments are all known sufficiently accurately [ |T2[ that the experimental errors 
play no role in our subsequent analysis.) We stress that these expressions for G%^ are exact 
consequences of QCD, under the assumption of charge symmetry. Equation (||) provides 
a particularly favorable case for the determination of with minimal dependence on 
the valence-quark ratio. 
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Figure 2. The consistency relation between Up/u-^+ and Un/u^o which must be satisfied 
within QCD. The part of the straight fine which is dashed corresponds to G%i{0) < 0, while 
the sohd part of the fine has G\.j{0) > 0. The standard quark model assumption of intrinsic 
quark moments independent of their environment is indicated by the filled square at (1, 1). 
The diamonds on the constraint line mark 0.1 fi^ changes in 6*^(0) when ^R^ = 0.55 [ 



T3[| . The lattice QCD prediction (after an appropriate chiral extrapolation, discussed in 



following sections) is illustrated by the filled circle with standard errors indicated. 



If one considers the quark model suggestions of Un/u^o = 1 and ^R^ = 0.65 in (P), 
one finds = —0.81 fiN, a significant departure from the experimental indication of 
positive values. 

Equating (J^) and (||) provides a hnear relationship between Up/uj2+ and Un/u^o which 
must be satisfied within QCD. Figure ^ displays this relationship by the dashed and solid 
line, the latter corresponding to values for which G%j{0) > when ^R'^ is in the anticipated 
range < ^i?^ < 1. Since the line does not pass through the point (1, 1) corresponding to 
the simple quark model assumption of universality, the experimentally measured baryon 
moments are signaling that there must be an environment effect exceeding 12% in both 
ratios or approaching 20% or more in at least one of the ratios. Moreover, a positive value 
for G%j{0) requires an environment sensitivity exceeding 70% in the Un/u^o ratio. Hence 
the experimental indication that G%j{0) > challenges the intrinsic magnetic moment 
concept which is fundamental to the constituent quark model. 
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3. CHIRAL CORRECTIONS 

One of the major challenges at present in connecting lattice calculations of hadronic 
properties with the physical world is that computational limitations restrict the accessible 
quark masses to values an order of magnitude larger than the physical values. At such 
large masses one is far from the region where chiral perturbation theory is applicable. 
Yet one knows that for current quark masses near zero there is important non-analytic 
structure (as a function of the quark mass) which must be treated correctly if we are 
to compare with physical hadron properties. Our present analysis of the strangeness 
magnetic form factor has been made possible by a recent breakthrough in the treatment 
of these chiral corrections for the nucleon magnetic moments [ |T1|, 0. In particular, a 



study of the dependence of the nucleon magnetic moments on the input current quark 
mass, within a chiral quark model which was fitted to existing lattice data, suggested a 
model independent method for extrapolating baryon magnetic moments which satisfied 
the chiral constraints imposed by QCD. We briefiy summarize the main results of that 
analysis: 



a series expansion of /ip(„) in powers of is not a valid approximation for larger 
than the physical mass. 



• on the other hand, the behavior of the model, after adjustments to fit the lattice 
data at large m,r is well determined by the simple Pade approximant: 

/^p(n) = ^ ; k ■, (9) 

1 — ^"^TT + cm^ 

as illustrated in Fig. ^. 

• Eq.(^ not only builds in the usual magnetic moment of a Dirac particle at mod- 
erately large but has the correct leading non-analytic (LNA) behavior of chiral 
perturbation theory 

= yUo + X^^TT, (10) 

with X ^ model independent constant. 

• fixing X at the value given by chiral perturbation theory and adjusting /io and c 
to fit the lattice data yielded values of /ip and fin of 2.85 ± 0.22 /xtv and —1.96 ± 
0.16 /xtv, respectively, at the physical pion mass. These are in remarkably good 
agreement with the experimental values - certainly much closer than the usual 
linear extrapolations in niq. 

Clearly it is vital to extend the lattice calculations of the baryon magnetic moments 
to lower values of than the 600 MeV used in the study just outlined. It is also vital 
to include dynamical quarks. Nevertheless, the apparent success of the extrapolation 
procedure gives us strong encouragement to investigate the same approach for resolving 
the strange quark contribution to the proton magnetic moment. 
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Figure 3. Extrapolation of lattice QCD magnetic moments (• Ref. [ P] and ■ Ref. [[T^) 
for the proton (upper) and neutron (lower) to the chiral limit. The solid curves illustrate 
a two parameter fit of Eq. (P) to the simulation data in which the one-loop corrected 
chiral coefficient of m.,r is taken from xPT. The dashed curves illustrate the leading non- 
analytic contribution of Eq. (plOl) . The experimentally measured moments are indicated 
by asterisks. 



Our focus here is to extrapolate existing lattice QCD estimates of valance quark contri- 
butions to baryon moments (up, Un, Me+, w=o) to the physical mass regime. The leading 
non-analytic pieces of these various contributions to the baryon magnetic moments are 
proportional to 

— 7^^^ = xm^, (11) 

where the pion decay constant /tt = 93 MeV. Isolating a particular quark flavour con- 
tribution only requires setting the electric charge of all other quark flavours to zero. 
Separation of the contributions of the LNA terms into valance and sea-quark sectors has 
been resolved in Ref. [[1^. The coefficients (3 and x for various quark sectors and baryons 
are summarized in Table |l|. In removing the u sea-quark contribution from the w-quark 
sector, one eliminates the contributions of the u sea-quark to the LNA term. Since the 
2;-component of angular momentum in the pion-baryon system is positive, the u contri- 
bution to the moment is negative. Hence its removal suppresses the non-analyticity in 
the neutron moment and enhances the non-analyticity in p, and S"*" where the pion 
cloud acts to enhance the u-quark sector contribution to these baryon moments. 

For the calculation of G\^{0), the LNA chiral behavior of m„ and is crucial, as 
the corresponding coefficients are of opposite sign. As a result, the ratio of these u- 
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Table 1 

Coefficients (3 and x of O and (11) for the leading nonanalytic (LNA) contribution of the 



M-quark(s) to the magnetic moment of various baryons. We express the LNA coefficients 
X in terms of the usual SU(6) constants F and D, as well as giving numerical values for x 
and j3 (which are equivalent through Eq.(^)). The w-quark charge has been normalized 
to 1. 



Baryon 




M-fiavor sector 


M-sea-quark loop 


M- valence sector 


n 


P 


{F + Df 


{9F^ -QFD + 5D^)/3 


2{-3 F^ + QFD - D^)/3 






1.020 


0.612 


0.408 




X 


4.41 


2.65 


1.76 






-{F-Df 


+ {F-Df 


-2{F-Df 




P 


-0.0441 


0.0441 


-0.0882 




X 


-0.191 


0.191 


-0.381 


P 


p 


-{F + Df 


{9F'^ -QFD + 5D'^)/3 


-4(3F2 + 2L)2)/3 




p 


-1.020 


0.612 


-1.632 




X 


-4.41 


2.65 


-7.06 


S+ 




-2(3F2 + D2)/3 


2(3F2 + D2)/3 


-4(3F2 + Z}2)/3 






-0.568 


0.568 


-1.136 




X 


-2.46 


2.46 


-4.91 



quark contributions, extrapolated with the correct chiral behavior, will be completely 
different from the linearly extrapolated ratio. That the signs of these chiral coefficients 
are opposite may be traced back to the charge of the predominant pion cloud associated 
with these baryons. The chiral behavior in the neutron case is dominated by the transition 
77. — > p7r~ ^ n, where there is a contribution from a u loop. On the other hand, for the 
S° the LNA u-quark contribution is dominated by the process H° 'E~n~^ S^. In this 
case the virtual 7r+ involves a valence w-quark - and hence the sign change. 

In order to extrapolate the lattice data for valence quarks to the physical pion mass 
we use the same approach which worked so well for the neutron and proton moments in 
Ref. [ pi |T6|. That is, we fit the lattice data for the valence quark of flavor q in baryon 



B with the form: 

Here x a^i'fi the modified LNA chiral coefficients for the valence quarks given in Table |l| 
(last column) and and c are fitting parameters. 

Figures ^ and |^ illustrate the extrapolations of the lattice data. The results of these ffis 
shown at = 140 MeV are the results of the extrapolation procedure, including ffiting 
errors. Naive linear extrapolations are also shown to help emphasize the dominant role of 
chiral symmetry in the extrapolation process. 

Clearly the effect of the chiral corrections on the extrapolated values of the key magnetic 
moment ratio m„,/mho at = 140 MeV is dramatic. Instead of the value 0.72 ± 0.46 
obtained by linear extrapolation one now finds 1.51±0.37. It is vital that this ratio is much 
more consistent with the range of values found necessary to yield a positive value of G^, 
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Figure 4. Extrapolation of the valence w-quark magnetic moment contributions to the 
neutron (solid symbols) and 5" (open symbols). Naive linear extrapolations (dashed 
lines) are contrasted with the chiral extrapolations incorporating the leading nonanalytic 
behavior. The lattice data is taken from Ref. [ |^. 



as illustrated in Fig. |^ - values which seemed quite unreasonable in the constituent quark 
model. By comparison, the ratio Up/uY.+ at = 140 MeV is quite stable, changing 
from 1.14 ± 0.08 in the case of linear extrapolation to 1.14 ± 0.11 when the correct 
chiral extrapolation is used. The lattice point illustrated in Fig. |^ is consistent with 
the constraint obtained from charge symmetry. This suggests that any charge symmetry 
breaking in the experimentally measured moments is small. 

The effect of the chiral extrapolation on the ratio m„/m=o is the most dramatic, changing 
it from 0.72 ± 0.46 to 1.51 ± 0.37. Thus Eq. (^, which seemed certain to guarantee 
< 0, no longer does. Indeed, using the recent estimate [ of ^i?^ = 0.55, we obtain 
= —0.16 ± 0.18/i7v While this is still negative it does permit a small positive value 
within the error. 
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